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Abstract. We discuss the stability of the photometric redshift estimate obtained with the SED fitting method 
with respect to the choice of the galaxy templates. Within the observational uncertainty and photometric errors, 
we find satisfactory agreement among different sets of theoretical and empirical templates using the Hubble Deep 
Field North as a target galaxy sample. Our results suggest that, especially at high redshift, the description of 
the physical processes of photon absorption in the interstellar and intergalactic medium plays a dominant role in 
the redshift estimate. The specific choice of the template set, as long as this includes both normal and starburst 
galaxies, is in comparison a minor issue. 

Key words. Galaxies: distances and redshifts - Galaxies: evolution - intergalactic medium - dust, extinction - 
Methods: data analysis - Techniques: photometric 



1. Introduction 

Physical properties of high-redshift galaxies have been ex- 
tensively explored in recent years thanks to a number of 
deep spectroscopic surveys carried out with the largest 
ground optical telescopes (e.g. Lilly et al. 1995; Cowie et 
al. 1996; Ellis et al. 1996). These surveys have enabled to 
trace well galaxy evolution up to z ~ 1.5. 

Beyond this redshift limit, however, ground-based 
spectroscopy is forcedly confined to the brightest objects, 
a few magnitudes above the limit achieved by deep pho- 
tometry For the time being, spectra of galaxies fainter 
than / ~ 24 AB mag are still difficult to obtain, and one 
must therefore rely on multicolor observations coupled to 
photometric redshift techniques to estimate redshifts for 
the majority of galaxies in this magnitude range. 

Photometric redshift techniques have already been ap- 
plied by many authors to the ultra-deep images in the 
Hubble Deep Field North (HDFN, Williams et al. 1996), 
addressing in fine detail the study of galaxy luminos- 
ity function, star formation rate (SFR), UV cosmologi- 
cal density, and clustering properties (Madau et al. 1996; 
Connolly et al. 1997; Sawicki et al. 1997; Franceschini et al. 
1998; Pascarelle et al. 1998; Connolly et al. 1998; Miralles 
& Pello 1998; Arnouts et al. 1999). 

Send offprint requests to: M. Massarotti 



Two main approaches have been used in the literature 
to measure photometric redshifts: the so-called training- 
set method (e.g. Connolly et al. 1997; Wang et al. 1999), 
and the spectral energy distribution (SED) fitting method 
(e.g. Sawicki et al. 1997; Giallongo et al. 1998; Bemtez 
2000). 

In the training-set method the relation between col- 
ors (sometimes also magnitudes) and redshift is calibrated 
using a spectroscopic galaxy sample. The accuracy of this 
method depends on the number and completeness of the 
reference templates, and obviously increases when acquir- 
ing new spectroscopic data. It becomes less reliable at 
high redshifts because of incompleteness and increasing 
redshift uncertainty of the most distant reference galaxies 
(see Bemtez 2000 and Fernandez-Soto et al. 2000 for a 
full discussion of the method drawbacks). 

In the SED fitting method, on the other hand, photo- 
metric redshift is obtained by comparing observed galaxy 
fluxes, /° bs at the i-th photometric band, with a library of 
reference fluxes, /* empl (z, T), depending on redshift and a 
suitable set of parameters T, that account for galaxy mor- 
phological type, age, metallicity, dust reddening etc. In 
this way, for each galaxy one can perform a x 2 confidence 
test obtaining the values of z and T that minimize flux 
residuals between observations and reference templates. 
In addition to the best-z estimate, the method identifies 
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an "optimum template" , supplying information on galaxy 
age and morphology, although special caution must obvi- 
ously be taken in using this piece of information, because 
of possible numerical degeneracies in the parameter space. 

In this paper we will discuss in some detail the perfor- 
mance of the SED fitting method as a function of the tem- 
plate set chosen to reproduce galaxy colors. The primary 
aim is to assess the effects of different template proper- 
ties on the numerical stability and physical consistency of 
the method. In our opinion, this "sanity check" must be 
preliminary to any direct comparison of photometric and 
spectroscopic redshift determinations, as different tem- 
plate sets can give statistically similar results when ap- 
plied to the same spectroscopic subsample, but produce 
a different redshift distribution on the entire photometric 
catalog (Yee 1998). 

The plan of the paper is as follows. We will first recall 
in Sec. 2 some key steps in the method. The theoretical 
and empirical template libraries are introduced in Sec. 3, 
while Sec. 4 presents results obtained by comparing the 
different galaxy template sets. The role of physical pro- 
cesses in the interstellar and intergalactic medium is dis- 
cussed in Sec. 5, where we show their relevant impact on 
redshift estimate. Sec. 6 finally summarizes our conclu- 
sions. 

2. The SED fitting algorithm 

Before proceeding with our discussion, it might be useful 
to briefly remind how the SED fitting algorithm works. 

As a first step one has to convert galaxy observed mag- 
nitudes for each i-th photometric band into incoming ap- 
parent flux, /° bs (A). This is equivalent to reconstruct the 
SED of target galaxies at very low spectral resolution by 
sampling their luminosity at the effective wavelength of 
the photometric bands available. 

For each object, the sampled flux has then to be com- 
pared with the reference spectral libraries of template 
galaxies, /' empl (z, T), computing a \ 2 merit function of 
the fitting residuals such as 

z=l °i 

In the equation, N is the number of photometric bands, 
and <7i is the observational uncertainty in the i-th. band. 
The scale factor s is chosen in such a way as to minimize 
X 2 for each template: 

_ EL/f bs /r p (^)M 2 m 

We will search the galaxy redshift in the range z G 
[0.0, 5.0] with a step of 0.02. 

As in the current literature, a further constraint has 
been added to the fit imposing that the age of the best 
template cannot be older than the Hubble time at the 
selected redshift (in our calculations we assumed a H = 



50kms _1 Mpc -1 , qo = 0.0 cosmology with A = 0). It can 
be verified, however, that redshift estimates do not depend 
critically on the details of the cosmological model, and this 
constraint could even be dropped with no serious impact 
on the final output. 

3. Theoretical vs. empirical galaxy templates 

A key issue in the SED fitting method is the selection of 
the reference flux library. In particular, one has to decide 
whether it is more appropriate to use empirical or syn- 
thetic galaxy templates. 

The major advantage of the first choice is that the ob- 
served SEDs for a suitable set of local galaxies, spanning 
the whole range of Hubble morphological types, gives by 
definition a physically consistent picture of the real galax- 
ies, at least the nearby ones. On the other hand, there 
is no obvious argument supporting a straightforward ex- 
tension of the current galaxy properties to high-redshift 
objects. At earlier cosmological epochs, evolution could 
play a substantial role in changing both morphological 
and spectrophotometric properties of distant galaxies (see 
Buzzoni 1998 for a critical discussion on the subject). 

In this respect, synthesis models should be preferred, 
as they do try to take into account evolution. Of course, 
some differences exist among current theoretical codes, es- 
pecially in the SED predictions for the UV range (Chariot 
et al. 1996), and these differences should be carefully ex- 
amined when trying to match high-redshift observations 
to derive z from multicolor photometry. Optical (or even 
infrared) observations of high-redshift objects sample UV 
restframe emission. 

We should also recall that at high z the restframe wave- 
length range explored by a given photometric system be- 
comes narrower. For example, at z ~ 1 the whole HST 
WFPC2 photometric system (i.e. U300, £?450) ^606 phis the 
/814 band) is just probing the wavelength range shortward 
of the galaxy 4000 A Balmer break, while for z ~ 2 the 
spectral window shrinks even more to a scarce ±800 A in- 
terval around 1800 A in the galaxy restframe (cf. Fig. [j]). 
Note in addition that, at least for early-type galaxies, SED 
drops to nominal values shortward of the Balmer discon- 
tinuity, so that k corrections become increasingly positive 
with z. As we can see in Fig. [lj the effect can reach up 
to two orders of magnitude for elliptical galaxies, and this 
can result in a strong bias against their detection at high 
redshifts. 

Another problem to be addressed when comparing 
photometric observations of high-redshift galaxies with 
reference templates is the role of intervening interstellar 
(ISM) and intergalactic medium (IGM), which strongly 
modulates the observed galaxy SED. Interstellar dust at- 
tenuates the galaxy luminosity by an amount that greatly 
depends on the specific galaxy environment (Gordon et 
al. 1997; Bruzual et al. 1988; Kuchinski et al. 1998). The 
IGM, on the other hand, drastically absorbs galaxy emis- 
sion shortward of the Ly-a limit (the well known Ly- 
a forest effect, e.g. Madau 1995). Unfortunately, present 
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2000 4000 6000 8000 10 4 
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Fig. 1. Restframe portion of galaxy SED explored by the 
WFPC2 photometric system at z = and z — 2. As an 
illustrative example a template model for local elliptical 
galaxies is shown. 

knowledge of these physical processes does not yet per- 
mit any firm and univocal correction of the data, leaving 
room for a variety of operational approaches (Sawicki et 
al. 1997; Giallongo et al. 1998; Miralles & Pello 1998). 

In our study we will include three sets of theoreti- 
cal galaxy models, namely those provided by the code of 
Bruzual & Chariot (1993, hereafter BC), Buzzoni (1989, 
1995, 2000; BUZ) and Fioc and Rocca-Volmerange (1997; 
FRV), as well as the empirical template set of Coleman 
et al. (1980; CWW) extensively used by previous authors 
(Sawicki et al. 1997; Fernandez-Soto et al. 1999, hereafter 
FSLY; Benftez 2000; Thompson et al. 2000). 

3.1. Template libraries 

For each galaxy we are mainly interested in a measure 
of its redshift, but an age estimate and a morphological 
classification parameter are also a result of the SED fitting 
technique. In order to ensure an homogeneous sampling 
of the galaxy distinctive parameters, a selected grid of 
reference models for each of the three synthetic libraries 
has therefore been considered. 

According to the leading physical processes in the 
galaxy stellar populations we can identify two main evolu- 
tionary regimes that induce relevant differences in galaxy 
luminosity: 

i) "Early" galaxy evolution - In the first Gyr of life, 
most of the galaxy luminosity is provided by short-living 
stars of high mass [M > 2M Q , Buzzoni 2000). These are 
strong UV emitters and dominate galaxy luminosity at 
short wavelength. Given their short lifetime, these stars 
track actual star formation activity in the host galaxy 
through a measure of its integrated UV flux. Synthesis 
models all converge toward a unified starburst scenario 
for t < 1 Gyr, despite any later morphological evolution 
that could then differentiate among spirals, ellipticals and 
irregulars (Buzzoni 1998; Leitherer et al. 1999). 



ii) "Late" galaxy evolution - It comprises all the ob- 
jects older than 1 Gyr, displaying morphological fea- 
tures accounted for by the standard Hubble classification. 
Disregarding any specific detail of star formation and mor- 
phology, the stellar bulk in old galaxies is dominated by 
stars with M < 2M Q experiencing a full red-giant-branch 
(RGB) phase tipping with a degenerate-core Helium flash 
(Renzini & Buzzoni 1986). 

For each of our synthesis codes, we first set up an "old- 
galaxy" template library, tuning up reference models for 
present-day galaxies along the whole Hubble sequence, 
and consistently tracking back-in-time evolution up to 
t = 1 Gyr. For ages less than 1 Gyr all our synthesis 
codes have been complemented by the starburst models at 
t = 50, 100, 500, and 800 Myr from Leitherer et al. (1999), 
better suitable to describe high-mass stellar evolution. 

Note that having used Leitherer et al. (1999) starburst 
models to describe young galaxies does not influence our 
results, as all the codes do actually converge, for early 
ages, to the redshift results provided by the Leitherer et 
al. (1999) models. 

An important issue to optimize computing resources 
and prevent redundant oversampling in the parameter 
space is the age step in the grid of template models. 

Based on the theory of simple stellar populations (cf. 
e.g. Buzzoni 1989) we know that luminosity evolution pro- 
ceeds in this case like L cx t~@ , with j3 varying in the range 
between unity at UV wavelength and 0.7 in the infrared 
(K band). Even accounting for a smoother evolution, e.g. 
for a continuous SFR, a A t/t ~ 0.3 age sampling is suf- 
ficient to confidently appreciate any 0.1 mag change in 
apparent colors for high-z spirals and ellipticals. 

The adopted BC templates consist of models with four 
different SFR histories to match today's E/S0, Sa/Sb, 
Sc/Sd, and Im Hubble types. From the BUZ and FRV 
codes, we selected the original templates for E, Sb, Sc, 
Im morphological types. The BC models assume a Scalo 
IMF, while the FRV models assume a Rana & Basu (1992) 
mass distribution, and the BUZ models adopt a standard 
Salpeter IMF. For all the models we sampled the differ- 
ent Hubble types at t = 1.0, 1.4, 2.0, 4.0, 8.0, 10.0, 12.0 and 
15.0 Gyr. 

Table |l| gives a complete summary of the relevant dis- 
tinctive properties of each theoretical library, reporting 
also the expected colors for present-day galaxies of differ- 
ent Hubble type. 

The original CWW empirical template set consisted 
of four observed galaxy spectra, E, Sbc, Seal and Im re- 
spectively. Sawicki et al. (1997) and Bcmtez (2000) had 
already noticed that adding starburst templates (sensibly 
bluer than any Im SED) to CWW greatly alleviates the 
systematic discrepancies found in redshift distribution of 
target galaxies beyond z > 1.5. Similarly, the lack of star- 
burst templates could have biased FSLY results by induc- 
ing an artificial peak in the galaxy distribution around 
z ~ 1.7. In this work, following the strategy adopted by 
Sawicki et al. (1997), we extended the CWW empirical li- 
brary by adding two starburst models from Leitherer et al. 
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Table 1. Galaxy templates from the reference synthesis codes and from Coleman et al. (1980) 



Template Age U - V B -V V - K SFR IMF<*' Remarks 

(* = 0) 



Bruzual & Chariot (1993): 



E/SO 


12 Gyr 


1.58 


1.00 


3.22 


oc exp(-t /t), t 


= 1 Gyr 


Sea 


Z = 




Sab 


12 Gyr 


1.09 


0.81 


2.95 


oc exp(-t/r), r 


= 4 Gyr 


Sea 


Z 


z e 


Scd 


12 Gyr 


0.78 


0.66 


2.74 


oc exp(-t/r), r ■ 


= 15 Gyr 


Sea 


/ 


z e 


Im 


1.4 Gyr 


0.13 


0.28 


1.89 


const. 




Sea 


/ = 


z 



Buzzoni (2000): 

E 15 Gyr 1.48 0.94 3.27 single burst Sal two-zone model, 

Lbol (halo:bulge) = (0.15:0.85) 
Z(halo + bulge) = Z e 

Sb 15 Gyr 0.58 0.63 2.79 oc t~ - 5 (disk) Sal three-zone model, 

Lbol (halo:bulge:disk) = (0.05:0.30:0.65) 
Z(halo + bulge) = Z e , Z(disk) = Z Q /3 

Sc 15 Gyr 0.36 0.53 2.62 oc t~° 1 (disk) Sal three-zone model, 

Lbol (halo:bulge:disk) = (0.03:0.15:0.82) 
]?(halo + bulge) = Z e , Z(disk) = Z Q /3 

Im 15 Gyr 0.16 0.43 2.39 oc f+° 8 Sal Z = Z Q /3 



Fioc & Rocca Volmerange (1997): 



E 


15 Gyr 


1.40 


0.91 


3.19 


Schmidt law 


R&B 


Schmidt parameters (n, 
Z = Z Q 




= (1, 0.2 Gyr), 


Sb 


12 Gyr 


0.92 


0.73 


3.14 


Schmidt law 


R&B 


Schmidt parameters (n, 

z = z e 




= (1, 2.8 Gyr), 


Sc 


12 Gyr 


0.59 


0.58 


2.94 


Schmidt law 


R&B 


Schmidt parameters (n, 

z = z e 




= (1, 10.0 Gyr) 


Im 


4 Gyr 


0.20 


0.36 


2.26 


Schmidt law 


R&B 


Schmidt parameters (n, 

z = z e 




= (1, 20.0 Gyr) 



Coleman et al. (1980): 

E - 1.49 0.93 - - - M31 model 

Sbc - 0.59 0.58 - - - - 

Scd - 0.37 0.49 

Im - -0.06 0.32 - - - - 



Leitherer et al. (1999): 



Starburst #1 


800 Myr 


-0.27 


0.13 


1.62 


const. 


Sal 


Z = Z e 


Starburst #2 


500 Myr 


-0.37 


0.09 


1.62 


const. 


Sal 


z = Z Q 


Starburst #3 


100 Myr 


-0.66 


0.02 


1.70 


const. 


Sal 


Z = Z Q 


Starburst #4 


50 Myr 


-0.76 


0.00 


1.72 


const. 


Sal 


Z = Z Q 



<*' Sal = Salpeter (1955), Sea = Scalo (1986), R&B = Rana & Basu (1992). 



(1999) at t = 50 and 500 Myr. This enlarged dataset will 
hereafter be referred to as the CWW-extended (CWWE) 
set. 



3.2. ISM and IGM corrections 

Fluxes provided by synthetic SED models have to be cor- 
rected for ISM and IGM absorption effects. Both are wave- 
length dependent, the latter is also redshift dependent, 
and they can dramatically change the galaxy spectrum in 
the UV restframe. 

Interstellar dust attenuation depends on the color ex- 
cess parameter E(B — V), scaling with dust column den- 
sity, as well as on the differential ISM opacity fc(A), which 



is related to the physical properties of dust grains. The 
classical relationship for intrinsic and emerging flux is: 

F obs (A)=^ (A)10-°- 4 ^- l/ )^). (3) 

In our work we adopted the Calzetti (1999) dust at- 
tenuation law fc(A) (see Fig. ||, panel a) including grain 
absorption and scattering. The E(B — V) parameter can 
take the values 0.0, 0.05, 0.1, 0.2, 0.3. By definition, the 
empirical CWW spectra already include internal redden- 
ing. 

We describe the IGM absorption as a function of red- 
shift following Madau (1995) (Fig. §, panel b). 

To give an illustrative example of the combined action 
of dust and IGM, in Fig. ||the BUZ Sb template has been 
reddened by different amounts of the color excess E(B—V) 
up to 0.3 magnitudes (upper panel). The effect of the IGM 
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2000 4000 6000 8000 
Wavelength [A] 

Fig. 2. Panel a: the adopted dust attenuation law after 
Calzetti (1999). Panel b: the IGM attenuation in a mag- 
nitude scale at redshift 2.5, 3.5 and 4.5, as labelled, ac- 
cording to Madau (1995). 



on the same spectrum is displayed on the contrary in the 
lower panel for the two cases of redshift z = 1.5 and z = 3. 

The galaxy spectral break induced by the cosmological 
Ly-a forest turns out to be the most effective feature to 
estimate redshifts beyond z ~ 2. The so-called "dropout" 
technique has been widely used to select high-redshift ob- 
jects in current deep surveys (Steidel et al. 1995; Madau et 
al. 1996). Notice though that for faint objects the strong 
signature due to the IGM external feature can be confused 
with the 4000 A Balmer break (which, on the contrary is 
intrinsic to galaxy SED) resulting in a biased redshift es- 
timate. 

4. Matching the FSLY catalog 

The HDFN photometric catalog produced by FSLY will 
be used as a target dataset for our analysis. It has been 
widely adopted in the recent literature (see, for instance, 
Arnouts et al. 1999; FSLY; Bemtez 2000) and consists 
of 1067 objects observed with HST WFPC2 in the four 
photometric bands L/300, -B450, V606> an d isi4) an d with 
the IRIM camera at KPNO in the infrared range (Johnson 
J,H, and K, Dickinson et al. (2000). The FSLY catalog 
contains fluxes, /° bs , and errors, cr^, for each galaxy in 




1000 10 4 

Wavelength [A] 

Fig. 3. The effect of internal reddening and IGM on the 
Sb reference template from Buzzoni (2000). Dust attenu- 
ation for E(B — V) up to 0.3 mag, as labelled, is shown in 
the upper panel, while the expected break induced by the 
Ly-a forest at z = 1.5 and z — 3 is shown in the lower 
panel. For reference, the HST photometric system and the 
Johnson JHK bands are displayed at the bottom. 

all photometric bands (i — 1, ...7). These quantities are 
those used in eqs. (|l|)-(||). 

The availability of IR observations is of special rele- 
vance for our analysis because, even for galaxies at z ~ 4, 
it allows us to sample restframe emission at optical wave- 
lengths, longward of the Balmer break. This greatly im- 
proves our chances to correctly interpret galaxy colors and 
obtain a fair redshift estimate for the most distant ob- 
jects. As we can see in Fig. [j], within the photometric 
errors, both elliptical and blue star-forming galaxies at 
z ~ 2 show basically the same SED within the wavelength 
range covered by the HST photometric system. Only IR 
data can allow us to break color degeneracy. 

For our analysis we selected from the original FSLY 
dataset only those objects with good photometry, that is 
with S/N > 5 in at least two bands, in order to rely confi- 
dently on at least one color. Out of a total of 1067 galaxies, 
we were left with 1041 objects in our subsample. As a re- 
sult of our selection we will be biased against the most 
distant objects beyond z > 5. For these galaxies the Ly-a 
break enters the Vqoq band, making them detectable only 
in the Igu band (IR observations are not deep enough, in 
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Table 2. Statistical uncertainty in redshift estimate due 
to the photometric errors for each reference library (ten 
bootstrap simulations). 





Az 


z < 1.5 


z > 1.5 


Az| > 0.5 


|Az| > 1.0 


BUZ 


-0.01 


0.16 


0.24 


11% 


7% 


BC 


-0.01 


0.16 


0.23 


12% 


8% 


FRV 


-0.01 


0.17 


0.21 


11% 


7% 


CWWE 


0.00 


0.15 


0.26 


12% 


8% 



any case, to grant a good S/N detection for these distant 
objects). 

The 4000 A Balmer break being one of the most use- 
ful spectral feature for photometric redshift determina- 
tion, we will single out a "low- redshift" subsample, in- 
cluding those galaxies for which the break is still within 
the HST photometric bands (this happens for z < 1.5), 
and a "high-redshift" subsample (z > 1.5) for which the 
4000 A break enters the IR photometric bands. 

4.1. Internal uncertainty of reference libraries 

It is of special relevance for our aim to preliminarily as- 
sess the effect of the observational uncertainties on the 
accuracy of redshift determination with the SED fitting 
method. 

We used a simple bootstrap procedure to generate ten 
copies of the FSLY catalog, obtained by adding to each 
galaxy flux a correction randomly extracted from the gaus- 
sian distribution of its photometric error. 

For each of the reference template libraries we then 
compared, for each object,, the redshift obtained from 
the measured photometry (z p hot) and that from the boot- 
strapped galaxy photometry (z s i m ). 

We will eventually have: 



^ ~ ^ ^ (Zphot ^-sim) 



iVsi 



and 



^ ^ [-Zphot -^sim] 
Nairn 



(4) 



(5) 



where N s i nl = 10 x 1041 is the number of simulated galax- 
ies. 

Moreover, we will check for the fraction of "discrepant" 
(l-Zphot - Zsim| > 0.5) and "catastrophic" (|z p hot - %im| > 
1.0) outliers in the low- and high-redshift galaxy subsam- 
ples. 

Table || gives a summary of the bootstrap results ob- 
tained for each theoretical synthesis code as well as for the 
CWWE empirical templates. Figure ^ shows in addition 
the distribution of the redshift residuals. 



Since photometric errors increase at fainter magni- 
tudes, redshift uncertainty is more pronounced for the 
most distant objects: u z ~ 0.16 at z < 1.5, while a z ~ 0.25 
for the z > 1.5 subsample. The catastrophic outliers are 
always less than 10% of the global sample, and there are 
no systematic effects, as expected. 

In addition to the photometric errors, also physical de- 
generacies in the z vs. color space domain of the models 
could themselves result in a source of numerical instability 
for the redshift estimates. In this regard, different model 
sets could in principle provide also different "bias pat- 
terns" in the redshift distribution. 

Quite interestingly, Table § shows that this is not the 
case for the four template sets accounted here. This might 
indicate that all reference libraries contain the same kind 
of degeneracy. In other words the galaxy colors are inter- 
preted in the same way (from a statistical point of view) 
by the different template sets. 

In Figure || the distribution of redshift residuals for 
catastrofic outliers (| Az| > 1.0) has a well defined pattern, 
and it is actually possible to identify two main sources for 
this pattern. The first is the color degeneracy between red- 
shift 2^0 and 1.3 < z < 2.1, that, given the photometric 
bands available in FSLY, shows up in all the templates. 
In the redshift interval 1.3 < z < 2.1, the entire Lyman 
series has not yet entered the L^oo band (cf. Fig. ||), while 
the Balmer break is located between the Isu and J bands. 
The uncertainty in the photometric observations can eas- 
ily cause a featureless objects to move within the ranges 
z ~ and 1.3 < z < 2.1. The second major source of the 
pattern visible in Figure ^ are the biased z estimates oc- 
curing when, due to large photometric errors, the Balmer 
break is mistaken for the Lyman decrement or viceversa. 
When the photometric errors in the IR bands are large 
enough to diluite/enhance a real/fictitious Balmer decre- 
ment, the objects affected are shuffled between the redshift 
ranges < z < 1 and z > 2.5. 

4.2. Comparison among reference libraries 

We can now compare photometric redshifts obtained using 
different reference template libraries. 

Given two independent redshift estimates, namely 
■Zphot, l and Zphot, 2 from different libraries, we can again 
define a simple statistical moment of the residual distri- 
bution in the form: 



-t— _ (-Zphot,! — z phot, 2) 



Az = J2 



A. 



while a variance can be computed as 

[■Zphot, 1 ~~ Zphot, 2] 2 



"2 = £- 



(6) 



(7) 



where N g — 1041 is the number of galaxies. 

Also in this case, we should check for the fraction 
of discrepant (|z p hot,i — z p hot.2| > 0.5) and catastrophic 
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2 4 2 4 



Zphot (FKV) z phot (CWWE) 

Fig. 4. Comparison of the statistical uncertainty in the redshift estimate due to photometric errors for the BC, BUZ 
and FRV reference libraries, as well as for the CWWE empirical templates. Residuals are from a bootstrap simulation 
of the FSLY galaxy catalog (1041 objects for each of ten simulated catalog samples). 





Fig. 5. Comparison of HDFN expected redshift distribution, for different reference libraries. The solid line is for 
Az = 0.0, while long-dashed strip is for Az = ±0.5. The vertical strip (solid lines) in panel d encloses the catastrophic 
outliers coming from the misinterpretation of the Balmer/Lyman break (see text for further details). 
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Table 3. Comparison of redshift distributions from the 
different reference libraries 





Az 


z < 1.5 


z > 1.5 


Az > 0.5 


|Az > 1.0 


BC vs. FRV 


0.02 


0.09 


0.12 


3% 


3% 


BC vs. CWWE 


0.03 


0.12 
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6% 


5% 


BUZ vs. BC 


0.00 
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4% 


3% 


BUZ vs. FRV 
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(l-Zphotji — %>hot,2| > 1-0) outliers in the low- and high- 
redshift galaxy subsamples. 

In addition, we should also test whether Az is consis- 
tent with zero or there is some systematic drift in data 
residuals. Mean standard deviation should be explored 
against possible trends as a function of galaxy redshifts. 

Figure | and Table | summarize the relevant results. 
The agreement between BC, BUZ and FRV (panel a-c in 
Fig. ||) seems quite good, with less than 5% of catastrophic 
outliers, and a z below 0.15 over the entire redshift range. 

The comparison with the CWWE template set gives 
slightly poorer values for a z (cf. panel d in Fig. ||), espe- 
cially for the most distant galaxies in the FSLY sample, 
but again a small fraction of discrepant redshift objects. 

As we can appreciate in Table |3[ there are only 
marginal systematic drifts in the Az residuals. 

It should be noted that in the synthetic libraries the 
"primeval" evolutionary scenario is described with the 
same templates (the four starburst models from Leitherer 
et al. 1999). On the other hand, even by replacing, for ex- 
ample, the "primeval" BC library with the young galaxies 
provided by the BC code itself, results of the comparison 
with other models are left almost unchanged. 

Comparing the results of Table |^ with those of Table 
||, it is clear that the whole SED fitting procedure is largely 
insensitive to any specific choice of the model reference li- 
brary. Rather than depending on the detailed galaxy SED, 
the redshift is mainly constrained by the recognition of the 
Balmer and/or the Ly-a break in the galaxy flux distri- 
bution. 

As a quite extreme but instructive example, in Fig. ^| 
we tried a fit of the galaxy subsample with fiducial red- 
shift z > 2 (259 objects, as estimated in every one of our 
reference libraries) using as a template set two Kurucz' 
(1992) model atmospheres for stars of T = 27 000 and 
18 000 K. These temperatures are roughly consistent with 
a B0-B2 spectral type like for stars of 10-20 Mq that 
dominate starburst galaxy luminosity in the early 10-100 
Myr of life (Leitherer et al. 1999). The redshift distribu- 
tion relying on this "minimal" template set well agrees 
with that obtained using, for example, the FRV model li- 
brary (<r 2 = 0.18 and no catastrophic outliers). At high 



Fig. 6. SED fitting of the FSLY galaxy subsample beyond 
z > 2. The data distribution by matching the FRV tem- 
plate set is compared with a "minimal" reference library 
consisting of two Kurucz' (1992) model atmospheres for 
stars of spectral type BO and B2. The relative a z results 
0.18 with no apparent drift in the data distribution. The 
solid line is for Az — 0.0, while long-dashed strip is for 
Az = ±0.5. Sec text for a full discussion. 



redshifts, therefore, the SED fitting method provides ro- 
bust estimates of galaxy redshifts, even when using an 
extremely poor library of templates. It is clear, however, 
that such a simplified procedure cannot be of any utility 
if we are interested in a more accurate analysis of high 
redshift data, such as to investigate age distribution or 
galaxy morphology at different cosmological times. 

Another interesting check of internal consistency of the 
SED fitting method deals with the age-redshift constraint. 
We discarded such constraint, by allowing the procedure 
to choose also galaxy models older than the Hubble time. 
Plotting the galaxy age information and the redshift in- 
formation, both obtained from the photometric redshift 
algorithm, will allow to check for their reciprocal consis- 
tency. 

For this exercise we grouped HDFN galaxies such as 
to guarantee at least 100 objects per redshift bin. Fig. [?] 
reports our results according to the BC, BUZ and FRV 
codes. 

Over the whole sample, the percentages of the fitting 
templates older than the actual Hubble time {[H 0l q ] = 
[50,0.0]) with BC, BUZ and FRV codes are 3.9 %, 8.5 
% and 5.7 % respectively. For each theoretical code the 
empirical redshift-age relation is virtually compatible with 
any cosmological model with Hq = 50 and qo < 0.5. 



4.3. Catastrophic outliers 

Although wide general agreement exists among template 
libraries in predicting galaxy redshift, it is worth investi- 
gating in some detail the few striking cases of catastrophic 
discrepancies (i.e. |Az| > 1.0). 
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Fig. 7. The expected age-redshift relation from the BC, 
BUZ and FRV best templates by imposing no cosmological 
constraint. HDFN galaxies have been grouped in order to 
guarantee at least 100 objects per redshift bin. Solid, grid, 
and diagonal-shaded histograms are the 50%, 75%, and 
90% envelopes of galaxy age distribution. In each plot, 
solid and dotted lines are the theoretical age-z relation 
for Hq = 50 and qo = 0.0 and 0.5, respectively. 



There are basically two groups of such outliers, as in 
the distribution of redshift residuals due to photomet- 
ric errors. The first group consists of objects with fea- 
tureless spectrum, that is with no significant break sig- 
nature (within the photometric errors). This is typically 
the case of the faintest blue objects in the HDFN sam- 
ple. Due to color degeneracy between redshift z ~ and 
1.3 < z < 2.1, these galaxies are shuffled within the 
two intervals by different libraries. The redshift interval 
1.3 < z < 2.1 is intrinsically critical for the SED fitting 
algorithm given the photometric bands available to us. As 
already noted above, while the entire Lyman series has 
not yet entered the U300 band (cf. Fig. ||), the Balmer 
break is located between the is 14 an d J bands. Therefore, 
redshift estimate entirely relies on the recognition of the 
smooth SED, and little changes in the reference templates 
can easily lead to diverging data interpretations. 
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Fig. 8. The r/300— -B450 apparent color as a function of red- 
shift for the CWW Im empirical template (solid line) and 
Leitherer et al. (1999) 50 Myr starburst model (dashed 
line). In absence of any starburst template, a fraction of 
star-forming galaxies at z ~ 2 (with f/300 — -B450 ~ 0.8) 
might be interpreted as local (z < 0.2) irregulars. 



The second class of catastrophic outliers comprises 
those objects for which the Balmer break is mistaken for 
the Lyman decrement or viceversa, according to the tem- 
plate library chosen, generating a series of biased z esti- 
mates beyond z > 2.5. This effect is especially evident 
when comparing CWWE vs. BC (cf. panel d in Fig. |J). 
A too scanty sampling of galaxy morphology, like in the 
CWWE templates, is the cause of the relative exacer- 
bation of the effect. For example galaxies with no deep 
Balmer break and relatively blue optical colors can be in- 
terpreted as Lyman break galaxies, in absence of alterna- 
tive intermediate age elliptical templates. 

More in general, a major source of catastrophic er- 
rors in the SED fitting procedure can result from the 
lack of starburst templates in the galaxy reference library. 
Standard Hubble types alone (even including Im irreg- 
ulars) are not "blue" enough to reproduce the colors of 
the most active star-forming objects seen at high red- 
shift. In Fig. H we show the C/300 — -B450 apparent color 
of the CWWE Im template (solid line) and that of the 
Leitherer et al. (1999) 50 Myr starburst model (dotted 
line) as a function of redshift. Lacking starburst templates 
can introduce a systematic shift in the redshift estimate 
(Az ~ 0.3) for those galaxies with (t/300 — -B450) ~ 0.8. 
Even worst, in some cases one might erroneously confuse 
UV-enhanced galaxies at z ~ 2 with local irregulars sam- 
pled in their optical SED. When comparing FSLY pho- 
tometric redshifts with spectroscopic redshifts both this 
effects are evident. 



5. Starbursts and the role of ISM and IGM 

Dust reddening and IGM absorption have a large impact 
on the detection of starburst -galaxies at high redshifts, as 
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these effects can effectively masque the original starburst 
spectrum. 

In a simple experiment using the CWWE library and 
imposing E(B — V) = in our fit, about 20% of galaxies in 
the FSLY catalog can suitably be matched by (reddening- 
free) starburst templates. When dust attenuation is ac- 
counted for, and E(B — V) is left to vary as a free fit- 
ting parameter, this fraction raises to about 55% of the 
whole galaxy sample. A confident evaluation of the frac- 
tion of UV-enhanced galaxies at the different distances 
would have a pervasive impact in a number of cosmologi- 
cal problems, like for instance the study of the cosmic SFR 
(Madau et al. 1996; Steidel et al. 1999; Thompson et al. 
2000). 

To explore the selective influence of interstellar dust 
on the SED fitting procedure, we considered the FSLY 
subsample with fiducial z > 1.5 (as estimated in every 
one of our reference libraries) and compared it with the 
CWWE empirical library (see Fig. ^). 

Out of 348 selected galaxies, 341 can be fitted by star- 
burst templates, Scd spirals or later morphological types. 

In a dust-free fitting scheme (i.e. by assuming E(B — 
V) = 0) we derive the following relative partition among 
fiducial morphological types: 

[others : Scd : Im : Starbursts] = [7 : 34 : 160 : 147] 

with about 12% objects poorly fitted, that is at least 3-cr 
off in each band according to eq. (jl]). Incidentally, if we 
exclude from the template library the starburst models we 
obtain: 

[others : Scd : Im] = [7 : 43 : 298] 

but at the cost of almost doubling the fraction of poorly 
fitted objects (~ 24% of the total sample). 

If we leave E(B — V) as a free fitting parameter (in the 
range 0.0 < E(B - V) < 0.3), the quality of the fit im- 
proves remarkably (only 5% of objects deviates by 3-cr or 
more), and an even larger number of starburst candidates 
is obtained: 

[others : Scd : Im : Starbursts] = [7 : 25 : 41 : 275]. 

A dust-free fitting scheme also drastically increases the 
number of catastrophic outliers (~ 10% on the entire cat- 
alog, cf. Fig. [h]), and at z > 1.5 the redshift of reddened 
galaxies is overestimated by Az ~ 0.2 since IGM absorp- 
tion, which is an increasing function of redshift, is forced 
to mimic color reddening. 

In Fig. [| the scatter of the HST data (blueshifted 
according the Zphot value and normalized according the 
scaling factor s of eq. |2|) around the CWWE templates is 
shown directly. As we can see, at A > 1216 A the best tem- 
plates are an accurate description of high redshift galaxy 
photometry. Notice however that at A < 1216 A the ob- 
served fluxes deviate from the best templates: the exper- 
imental data do contain the IGM absorption, while the 
best templates do not. 
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Fig. 9. Restframe composite SED from the HST data for 
the 341 HDFN galaxies with fiducial z > 1.5 compared 
with the CWWE reference templates. From top to bottom 
solid lines display the starburst template with t — 500 Myr 
and E(B - V) = 0.0,0.05,0.1,0.2,0.3, the t = 50 Myr 
starburst model with E(B - V) = 0.0,0.05,0.1,0.2,0.3, 
the Scd and the Im CWW spectra. 

The effect of IGM at high redshifts is explored in 
Fig. |l[ Redshifts in the FSLY catalo gue are determined 
using the BUZ reference library with and without IGM 
taken into account. A drift clearly appears beyond z > 1.5 
when the Ly-a break enters the U300 band. The shift is 
generated by the wavelength difference between the 912 
A Lyman break in the galaxy spectrum and the various 
steps of the Ly-a forest. The fraction of outliers amounts 
to about 12% on the entire catalog and to ~ 21% of the 
z > 1.5 subsample. 

6. Discussion and conclusions 

We investigated the stability of photometric redshift es- 
timates in the SED fitting method with respect to the 
galaxy template set chosen. We used as a target dataset 
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Fig. 10. The selective influence of galaxy internal redden- 
ing on the SED fitting performance. The HDFN redshift 
distribution from the BC template set is displayed with 
and without taking into account for ISM absorption. This 
is done by leaving E(B — V) = or as a free fitting pa- 
rameter in the range 0.0 < E(B — V) < 0.3. Long-dashed 
lines in the plot represent Az = ±0.5. 



O 



6 




r, 


1 1 V »1 

A, 5 ' 

/ / / 
/***'/' 


'V 

/ / — ^ 

/ / 
/ / - 
/ 

/ 


4 




7 / 

/V* 


/ A/ y 

/ J&$*!/ / 
/ «a»gP / / 
*5^9j! / / 

W s 




2 


~j 


/ 


X / 

/ 







_ / j. 


/ 

' , , 1 


1 1 1 1 1 1 


r 







2 

Z phot 


4 

with IGM 


6 



Fig. 11. Comparison of the HDFN redshift distribution 
obtained from the BUZ reference library with and with- 
out taking into account the IGM absorption. Long- dashed 
lines in the plot represent Az = ±0.5. The solid line shows 
the expected upper limit in the redshift drift due to the 
the misinterpretation of the intrinsic Lyman break with 
the Ly-a forest effect. 

the FSLY catalog in the HDFN and as libraries of galaxy 
colors those produced by the theoretical codes of Bruzual 
& Chariot (1993), Buzzoni (2000) and Fioc and Rocca- 
Volmerange (1997), as well as the empirical template set 
of Coleman et al. (1980). This way our model reference 
framework comprises (and enlarges) the sample of tem- 
plate libraries used by almost all the literature on the field. 

The existence of non negligible photometric errors in 
the data is by itself an important source of uncertainty in 
the redshift estimates. Comparing the FSLY results with 



a number of bootstrap simulation of the same galaxy sam- 
ple we estimate that the internal uncertainty in the SED 
fitting performance amounts to a z <~ 0.16 at z < 1.5 and 
a z ~ 0.25 at z > 1.5, with ~ 8% of catastrophic outliers. 
This result is almost independent from the adopted model 
library and can be read as a lower limit to the intrinsic ac- 
curacy of the fitting method given the photometric errors 
of the FSLY catalogue. 

The external redshift uncertainty, found by compar- 
ing results obtained using different theoretical models, is 
typically a z <~ 0.11 at z < 1.5 and a z ~ 0.13 at z > 1.5, 
with less than 4% of catastrophic outliers. Comparing the- 
oretical models and the CWWE empirical set we obtained 
a z ~ 0.13 at z < 1.5 and a z ~ 0.18 at z > 1.5, with less 
than 5% of catastrophic outliers. Our results show that, 
given the photometric errors in the data, the SED fitting 
procedure is to a large extent almost insensitive to any 
specific choice of the model reference library. 

The lack of starburst templates in the model library, on 
the other hand, introduces a systematic bias in the redshift 
estimate (Az ~ 0.3) and a sensible increase of catastrophic 
outliers as UV-enhanced galaxies at z <~ 2 can be confused 
with local irregulars sampled in their optical SED. 

Not including in the fitting procedure ISM dust red- 
dening or IGM opacity has a dramatic impact on the dis- 
tribution of redshifts measured. 

A dust-free fitting scheme drastically increases the 
number of catastrophic outliers (~ 10% on the entire cata- 
log) and at z > 1.5 the redshift of dust reddened galaxies is 
overestimated by Az ~ 0.2. When dust attenuation is ac- 
counted for, the fraction of galaxies matched by starburst 
templates raises from <~ 20% (by imposing E(B — V) = 0) 
to - 55% of the whole HDFN sample. 

In an IGM-transparent fitting scheme a drift in the 
redshift estimates clearly appears beyond z > 1.5 and the 
fraction of outliers amounts to ~ 21% of the high-redshift 
subsample. 

Our results suggest that, especially at high redshifts, 
the description of physical processes of photon absorption 
in the interstellar and in the intergalactic medium plays a 
key role in the SED fitting method for interpreting photo- 
metric data. Any specific choice of the templates, as long 
as these describe both normal and starburst galaxies, has 
on the contrary a minor effect on the performance of the 
method. 
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